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A computation of the Isgur-Wise Universal Function using a lattice implementation of the heavy quark effective theory is described. 



In this talk we review the lattice implementation 
of the heavy quark effective theory, describe its 
application to the computation of the universal heavy 
quark form factor, and show the first simulation 
results. 

Two other calculations of this decay form factor 
treating the heavy quarks as ordinary quarks but with 
a small hopping constant, instead of using the lattice 
heavy quark theory are also being presented at this 
conference. [1] 

The Isgur-Wise limit is the situation in which the 
mass of heavy quark is taken to be much larger than 
Aq CD , which is the characteristic energy of the 
strong interactions. [2] One separates the quark 
four-momentum into a fixed part proportional to the 
heavy mass and a residual dynamical part 



P = Mv + p 
v 2 = 1 



(1) 



The unit four-vector v is referred to as the classical 
velocity. In the infinite mass limit the quark 
propagator greatly simplifies: 
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The mass of the heavy quark drops out, and the 
Dirac structure of the propagator becomes simply a 
projector on positive energy states. This is, of 
course, the origin of the famous heavy quark spin- 
flavor symmetry. The heavy quark limit can be 
formulated as a field theory, which facilitates both 
non-perturbative analysis, the focus of this talk, and 
a systematic treatment of corrections in the inverse 
heavy quark mass. One constructs a field for a 
heavy quark of velocity v by factoring the rapid 



phase proportional to the large mass from the 
canonical field, 



/z <v) (x) - e iMv - x 1 + Y ' v y(x) 



(3) 



There is an independent field for each classical 
velocity. [3] The Lagrangian for h (v \ 



S£ (v) = h (v) (x) iv-D h M (x) 



(4) 



called the heavy quark effective theory, is the 
starting point of our analysis. 

To incorporate the heavy quark effective theory 
into a lattice calculation, we Wick rotate it into 
Euclidean space. [4] The classical velocity, which is 
only a parameter, is complex when expressed in 
terms of Euclidean four-momenta: 



v - (/v ,v) 



(5) 



v = \/ 1 + v 2 



The reduced propagator for the heavy quark satisfies 
the Euclidean space equation 



-ivd § 1 (x) = 8(x) 



with the boundary condition 

S {v \x) = 



x 4 < 



(6) 



(7) 



In the presence of gauge interactions, the derivative 
is replaced by a gauge covariant derivative. In order 
to most simply incorporate the boundary condition 
on the lattice, we represent the time derivative in the 



reduced Dirac equation by means of a forward 
difference, but use a symmetrical centered first 
difference for the space derivative. The lattice 
equation for the reduced propagator is: 

v [ U(x,x+i ) S M (x+t,y)- S (v \x,y) ] 

^ [£/(jc,jc+A) S M (jc+A,y) (8) 

- [/(x,x-p) 5 <0 (x-p,y) ] = 8(x,y) 



With the lattice heavy quark theory described 
above, we are in a position to address the 
paradigmatic heavy quark process, B — > D * I v, the 
weak decay of one meson containing a heavy quark 
into another. 

B -> D/V , D*/V 

(9) 

(Dv'| cy m (1 -iy 5 )b\Bv} ~ £(vv') 



This is solved by simple forward recursion; no 
iterative techniques are needed. 

Even for non-zero v, only trivial computation is 
needed to construct the heavy quark propagator. 
This makes it entirely feasible to compute three-point 
functions, and one may envision optimizing 
physical applications by a careful choice of heavy- 
light particle wave functions or replacing the simple 
first difference by an improved one which eliminates 
the leading lattice spacing corrections. 

The lattice heavy quark theory has some 
unfamiliar features which are worth noting. The 
symmetrical first difference gives rise to doublers. 
However, unlike the usual fermion situation, their 
contributions automatically vanish in the zero 
spacing limit, without the necessity of adding a 
Wilson-like term. The reason is that they always 
occur in association with gluon modes whose 
momenta are at the edge of the Brillouin zone, and 
those modes have energies on the order of the 
inverse lattice spacing. This is easily seen in 
perturbation theory. [5] To leading order in the 
inverse heavy quark mass, there are no closed heavy 
quark loops. 

Another feature of the lattice heavy quark theory 
is that for residual momenta oriented oppositely to 
the classical velocity, the propagator grows with 
Euclidean time. The reason is that the falloff is 
governed by the difference between the actual energy 
of a mode and that of the zero-residual-momentum 
mode, which is v-p/v Q in the heavy quark limit. 

When this negative, the falloff becomes growth. The 
gluon modes always compensate for this growth in 
simulations, but in perturbation theory, consistency 
requires an appropriate choice of contour in the 
Euclidean energy plane. [5,6] 



All the dynamical information is contained in the 
function known as the Isgur-Wise function, or the 
heavy quark universal form factor. It is a function 
only of the dot product of the classical momenta, and 
is normalized to unity for forward decay. 

= i do) 

The amplitude is contained in the reduced three-point 
function 

g'^Wwo) (id 

~ (Trs(y,x) T D S^\x,z) r,S iv '\z,y) T B ) 

The Isgur-Wise function may be extracted from 
the three -point correlation by dividing out the 
normalized single meson propagators. Better yet, we 
may take advantage of the normalization in the 
forward direction and use the three-point function at 
v D - v B as normalization. The expression from 

which we extract the universal form factor is then 

^(W) 2 = 

nm G ivy \x Q ,y ,z ) G (v '' v \x ,y ,z Q ) (12) 
,>,>, G (v ' v \x ,y ,z ) G^\x ,y ,z ) 

This form has the advantage that the overall 
normalization of the three-point function drops out 
completely. Only those renormalizations that are 
classical velocity-dependent need be included. 

We evaluated the three -point functions using the 
lattices made publicly available by Bernard and Soni. 
There were an ensemble of 16 lattices of size 
16 3 x 24 with lattice coupling P = 5.7. The light 
quarks were represented by Wilson quark 
propagators with hopping constant K = .164. The 



heavy quark propagators were evaluated for four 
values of each component of the initial and final 
quark classical velocity: 0, .25, .50, and .75. That is, 
there were 4 6 = 4096 combinations of v B and v n . Of 



D' 



course, the many that were lattice rotations of each 
other were averaged to improve statistics, but even 
so there were more than a hundred different values 
of v B -v D , lying in the interval [1,1.822]. 

To improve the overlap of the meson wave 
function with the ground state, we applied multiple 
single-link smearing steps. Five steps gave 
statistically cleaner signals than none, but ten steps 
was not an improvement over five. 

In order to project out the mesonic ground states 
from the three -point function, one either needs 
excellent wave functions or large Euclidean time 
separations between the slices on which the heavy 
meson is created, emits a weak current, and is 
detected. Unfortunately, the statistical noise in the 
simulation overwhelms the signal when those 
separations are as large as 4 units of Euclidean time. 
The largest separations that provide a signal are 
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(13) 



At 2 units of Euclidean time separation, which 
corresponds to a greater contamination with higher 
mass states than at 3 units, the computed form factor 
is flatter. This presumably reflects Bjorken's sum 
rule, which states that the sum of the squares of the 
Isgur-Wise functions to all final states is 1. [7] 

The results of the simulation together with a 
comparison to the data on B — > D * / v from both 
ARGUS and CLEO are shown in Fig. 1. [8] The 
agreement with the experimental data is rather 
spectacular. A simple quadratic fit to the simulation 
results gives for the slope at the origin 



.95 



(14) 



The strengths of the lattice heavy quark theory we 
have applied here are the same as those of the Isgur- 
Wise limit itself. By eliminating the heavy quark 
masses before undertaking a lattice analysis, only the 
residual four-momentum of the heavy quark is 
constrained by the inverse lattice spacing, and so 
there is in principle no restriction on v • v . 



ARGUS data 
□ CLEO data 




Four-vector v.v' 

Figure 1 . Lattice calculation of the Isgur-Wise form 
factor. Verticle bars are statistical errors for each 
simulation point. The solid curve is a quadratic fit 
to the simulation results. 



The calculation described here can be improved by 
a more optimal choice of wave function, by 
including the velocity-dependent continuum [9] and 
lattice [6] renormalizations, and by the use of 
improved lattices and propagators from which the 
leading finite lattice spacing errors have been 
removed. We hope to have results of such a 
calculation by the next lattice meeting. 

REFERENCES 

1. C. Bernard, Y. Shen. and A. Soni, Phys. Letts. B 317 (1993) 
164; R. Kenway (UKQCD Collaboration), this conference 

2. N. Isgur and M.B. Wise, Phys. Lett. B 232 (1989) 113; B 237 
(1990) 527: M.B. Voloshin and M.A. Shifman, Yad. Fiz. 45 
(1987) 463; Sov. J. Nucl. Phys. 47 (1988) 511 

3. H. Georgi, Phys. Letts. B 240 (1990) 447 

4. J.E. Mandula and M.C. Ogilvie, Phys. Rev. D (Rap. Comm.) 
45 (1992) R2183; Nucl. Phys. B (Proc. Suppl.) 26 (1992) 459 

5. J.E. Mandula and M.C. Ogilvie, Nucl. Phys. B (proc. Suppl.) 
30 (1993) 447 

6. U. Aglietti, HEP-PH/9304274 (1993) 

7. J. Bjorken, in La Thuile Rencontres (1990) 583 

8. H. Albrecht et al. (ARGUS Collaboration), Phys. Letts. B 219 
(1989) 121; G. Crawford et al. (CLEO Collaboration), preprint 
(1993) 

9. A. Falk et al., Nucl. Phys. B 343 (1990) 1 



